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Abstract 

The definition of the covariant space-time averaging scheme for the objects (tensors, 
geometric objects, etc.) on differentiable metric manifolds with a volume n-form, which 
has been proposed for the formulation of macroscopic gravity, is analyzed. An overview 
of the space-time averaging procedure in Minkowski spacetime is given and comparison 
between this averaging scheme and that adopted in macroscopic gravity is carried out 
throughout the paper. Some new results concerning the algebraic structure of the 
averaging operator are precisely formulated and proved, the main one being that the 
averaging bilocal operator is idempotent iff it is factorized into a bilocal product of 
a matrix-valued function on the manifold, taken at a point, by its inverse at another 
point. The previously proved existence theorems for the averaging and coordination 
bilocal operators are revisited with more detailed proofs of related results. A number of 
new results concerning the structure of the volume-preserving averaging operators and 
the class of proper coordinate systems are given. It is shown, in particular, that such 
operators are defined on an arbitrary n-dimensional differentiable metric manifold with 
a volume n-form up to the freedom of (n — 1) arbitrary functions of n arguments and 1 
arbitrary function of (n— 1) arguments. All the results given in this paper are also valid 
whenever appropriate for affine connection manifolds including (pseudo)-Riemannian 
manifolds. 
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1 Introduction 

Space-time averaging procedures play an important role in modern physics because of, first of 
all, their relevance to deriving classical macroscopic theories. A well-known example of such 
a procedure is the space-time scheme developed for averaging out the microscopic Lorentz 
electrodynamics to derive the macroscopic Maxwell electrodynamics (see, for example, [[j]]- 
PJ )H Another important physical argument for considering space-time averaging procedures 
is that they are relevant to modelling the process of physical measurement. It is the space- 
time averages of physical fields that are known JTD|] , [|T1]] to have direct observational status 
and physical meaning. 

The space-time averaging procedure of classical electrodynamics utilizes explicitly the 
flat character of the Minkowski space-time manifold and its formulation is essentially based 
on the existence of Cartesian coordinates. In this connection, the following questions are of 
interest and importance: whether it is possible (a) to encode the properties of the averaging 
procedure in a covariant manner suitable for differentiable manifolds not equipped with 
a metric and a connection in general, and (b) to formulate the corresponding covariant 
averaging procedure for the objects which can live on such manifolds and to clear up its 
geometric meaning. Having replied (a) and (b), the next question of primary importance is 
(c) how the averaging procedure may be made compatible with a metric and a connection 
when the differentiable manifold possesses those structures. 

It should be pointed out here that the above problem of generalizing the flat space-time 
procedure for curved manifolds goes far beyond being simply an academic problem. One 
of its most important areas of applicability is the general theory of relativity where the 
space-time is a 4-dimensional pseudo-Riemannian manifold. As is well-known, there is not 
yet a satisfactory derivation of the Maxwell equations in general relativity, apart from the 
covariantization procedure where the partial derivatives of the special theory of relativity are 
replaced by covariant ones (see, for example, Jl^]). A physically motivated and mathemati- 



cally correct derivation of the macroscopic Maxwell equations by averaging out the general 
relativistic microscopic Maxwell-Lorentz equations is still lacking and, furthermore, the very 
foundations of microscopic electrodynamics in general relativity are not well-established yet. 
Another problem in general relativity where the availability of a space-time averaging proce- 



dure is of primary importance is the so-called averaging problem (see Pl,[fL3|-[115 for a review 
and discussion). Its main motivation comes from cosmology where Einstein's equations are 
usually utilized with a hydrodynamic stress-energy tensor without any satisfactory proof of 
why the left-hand side of the equations (the field operator) keeps the same structure while the 
right-hand side has been changed, or averaged, from a real discrete matter distribution (stars, 
galaxies, etc.) to a continuous one. The task here is to carry out a space-time averaging of 
Einstein's equations in order to understand the structure of the averaged (macroscopic) field 
equations and apply them to deal with the overwhelming majority of cosmological problems. 
A solution for this problem is also desirable in order to provide a rigorous basis for con- 
structing continuous matter models from discrete ones (this construction is based in modern 
cosmology mainly on phenomenological grounds). 

The goal of this paper is twofold: (1) to give precise formulations and detailed proofs of 
some new results concerning the algebraic structure of the averaging and coordination oper- 



Alternative approaches for deriving macroscopic electrodynamics which apply other averaging proce- 
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ators and the properties of the space-time averages within the so-called macroscopic gravity 
P|)[1l4|)||I^|-||l9| (the basics of this averaging procedure has been developed in pO^pl]]); (2) 



to revisit the previously shown existence of the bilocal operators of macroscopic gravity 



I4[],[|T(| with more detailed proofs of related results. The paper is organized as follows. Sec- 



tion 2 gives an overview of the space-time averaging scheme used in electrodynamics. The 
space-time averages of macroscopic gravity are defined in Section 3. The next two sections 
are devoted to the algebra of the averaging operator with a discussion of its algebraic prop- 
erties and a formulation of an important theorem stating that the averaging bilocal operator 
is idempotent iff it is factorized into a bilocal product of a matrix-valued function on the 
manifold, taken at a point, by its inverse at another point. The differential properties of 
the averages and the coordination bivector are summarized in Section 6. In Section 7 the 
existence theorems for the averaging and coordination operators are formulated and proved. 
The last Section 8 is devoted to the definition and properties of the proper coordinate sys- 
tems, the coordinates in which the averaging operators take their simplest form. It is shown, 
in particular, that such operators are defined on an arbitrary n-dimensional differentiable 
metric manifold with a volume n-form up to the freedom of (n — 1) arbitrary functions of 
n arguments and one arbitrary function of (n — 1) arguments. A summary on the meaning 
and content of the space-time averaging procedure of macroscopic gravity is given in the 
Conclusion. 

Though given for n-dimensional differentiable metric manifolds with a volume n-form, 
they are valid whenever appropriate for affine connection manifolds, including (pseudo)- 
Riemannian manifolds. All the results in this paper also hold generally for n-dimensional 
differentiable manifolds with a volume n-form possessing neither a metric nor a connection. 



2 Space-time averages in Minkowski manifold 

The space-time averaging procedure applied in electrodynamics is explicitly based on the 
flat character of the Minkowski space-time £ and the existence of the Cartesian coordinates 
(t,x a ),a = 1,2,3. The space-time averages are defined as follows 0-0. 

Definition 1 The average value of a tensor field p^(x), x G S, over a compact space region 
S and a finite time interval At at a supporting point (t, x a ) G At x S is 

W^ x ^ £ = a7v5 f M f g i$<t + + ■ (!) 

Here V$ is the 3-volume of the region S, which is usually taken as a 3-sphere of radius R 
around the point x a at the instant of time t, 

V s = [ d 3 x' . 
Js 

Formula (|I|) defines the average value (pp(t,x a ))£ at a point (t,x a ) on the manifold 
£ . In order to obtain averaged tensor fields and define its derivatives, it is necessary to 
make additional assumptions concerning the averaging regions S and the intervals At. This 
assumptions are usually made only tacitly (see, however, a discussion in U), or they are 
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(pp(t,x a ))£, (ii) all the regions S and time intervals At are typical in some defined sense - 
they are usually required to be of the same shape and volume, V$ = const and At = const, 
and related to each other by shifting along the Cartesian coordinate lines. These properties 
are very easily arranged by Lie-dragging of a region S and an interval At, chosen around a 
point (if:, x a ) G £, along the congruences of the Cartesian coordinate linesQto get a "covering" 
of the manifold (or its connected part) with a region of the same shape and volume and a 
time interval of the same length around each point of £. Consequences of the properties (i) 
and (ii) are, first of all, the uniqueness of the definition of the average field (pp(t, x a ))z, and, 
secondly, the commutation formulae between the averaging and the partial derivatives 

j t (v%t,x a )) e = (^(t,x«))s , ^-M(t^ a ))s = (^- a PW,x a ))e • (2) 

The essence of this approach lies in the exploitation of the calculational advantages of the 
Cartesian coordinates on a flat manifold. The Cartesian coordinates play a central role both 
in defining the averages (|l|) and in obtaining properties, such as (@), which will allow the 
averaging out of the partial (ordinary) differential equations under interest. It should be 
stressed that by requiring (i) and (ii) the averages (Pp(t,x a ))£ become local functions of 
(t, x a ), i.e. 

( d d d d 



dx a dt dt dx a 



(p«(t,x a )) £ = 0, (3) 



and, therefore, they are provided with proper analytical properties. Furthermore, the func- 
tional dependence of the average ([]]) on the averaging region and interval, S and At, becomes 
simply a parametric dependence on the value of the volume Vg and the length At. One can, 
therefore, apply the standard differential and integral calculus to deal with the averages and 
averaged equations. The set of averaged components (p^(t,x a ))£ given by ([!]) in Cartesian 
coordinates are the components of a Lorentz tensor within the class of coordinate transfor- 
mations x a = ApX@ + a 13 with a constant shift a@ and a constant Lorentz transformation 
matrix A^. Now, if it is necessary to consider an average field in another coordinate sys- 
tem, the components of (pp(t,x a ))£ in that new system are found by applying the standard 
tensorial transformation law. 

There is still another property of the averages (pp(t,x a ))g widely used in macroscopic 
electrodynamics (but again, often only tacitly assumed); the idempotency of the averages 

((p«(t,x a )) £ ) £ = (p ( l(t,x a )) £ . (4) 

In order to prove this property we must calculate the average value ((pp(t, x a ))s )g of the 
average (pp(t,x a ))£ 

(W. *°))')' = *k L L (aW L L + ( ' + x " + x " + ^"^") «*>> • 

(5) 

Now, the expression ([|) leads to (^) under either of two additional assumptions: (Hi') the av- 
eraging region At' x S' is the same region At x S with the supporting point at (t', x' a ) G At x S 

2 A definition of Lie-dragging, or dragging, of a region of a manifold along a vector field (congruence), 
or in other words, mapping of a region into another along a vector field (congruence) can be found in any 
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and the average value (pp(t, x a ))g does not depend on the choice of a supporting point within 
a chosen region; (Hi") the averaging region At' x 5" is a different neighbouring region and 
the average value {pp(t,x a ))s remains the same if evaluated over any neighbouring regions 
containing a neighbourhood of x. Although they seem different the above assumptions are 
nevertheless essentially equivalent. Indeed, (Hi ') puts emphasis on the independence of the 
average value with respect to the choice of a supporting point from a set of all possible points 
in a fixed averaging region, while (Hi") puts emphasis on the independence of the average 
value with respect to the choice of an averaging region from a set of all possible regions 
defined by a neighbourhood of a fixed supporting point. Remembering (ii) that the averag- 
ing regions are typical, change of the supporting point can be considered as change of the 
averaging region, and vice versa. Both assumptions thereby encode the same fundamental 
property of the averages which is inherently related to the philosophy of averaging itself - 
an averaging region is considered as a point from the macroscopic point of view, and the 
change of either a supporting point or an averaging region does not affect the corresponding 
average value. The reason for this is that the differences caused by the above variations are 
negligibly small from a macroscopic point of viewf]. 

Thus, the well-known procedure of space-time averaging in classical electrodynamics 
presupposes the above specific conditions (i), (ii) and (Hi) to ensure reasonable analytical 
and tensorial properties of the averages ([!]). Clearly, the whole procedure relies crucially 
on the existence and properties of the exceptional coordinate system in a flat space-time, 
namely, Cartesian coordinates. These conditions, which seem more or less trivial at first 
sight, require further analysis to make clear their geometrical meaning and invariant content. 
This will allow a reasonable formulation of a space-time averaging procedure on general (not 
necessarily flat) manifolds. 



3 Definition of space-time averages in macroscopic grav- 
ity 

Let us remind the definition of the space-time averages adopted in macroscopic gravity 
|TJj],[]nj. This procedure is a generalization of the space-time averaging procedure adopted 
in electrodynamics (see Section 2) and it is also based on the concept of Lie-dragging of 
averaging regions, which makes it valid for any differentiable manifold. 



Definition 2 Chosen a compact region E C Ai in an n- dimensional differentiable metric 
manifold (Ai, g a p) with a volume n-form and a supporting point x G X to which the aver- 
age value will be prescribed, the average value of an object (tensor, geometric object, etc.) 

3 It should be noted here that in classical hydrodynamics, as distinct from macroscopic electrodynamics, 
discussion on the definition and properties of the averages continues for more than one hundred years. A 
definition of an average (either over space, time, ensemble, or a combination of such) in hydrodynamics and 
its properties are vital elements of the theory itself for it is clearly understood that the form of the equations 
depends on the definition and properties of the average. The definition (|l|) under conditions (i) and (ii) 
with the properties (^|), (|^) and (^), which are part of the Reynolds conditions in hydrodynamics, is known 
to result in the fundamental equations of hydrodynamics describing the dynamics of turbulence. If one of 
the Reynolds conditions is absent one must get different equations. For a discussion on averages and their 
properties in hydrodynamics, see, for example, p4j and references therein. 
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Pp(x), x E M. , over a region £ at the supporting point x G S is defined as 

= ^ jf p|(x, x')^g~'d n x' = <p«> , (6) 
where Vs is the volume of the region S ; 

= / e v 73 ^ . (7) 

Here the integration is carried out over all points x' G E, g' = det(g a p(x')) (since the 
primary interest is in space-time manifolds, it is assumed that this determinant is negative, 
otherwise the negative sign in \J— g must be dropped) and the bold face object pp(x,x') in 
the integrand of (|B|) is a bilocal extension of the object Pp(x), 

p^(x,x f )=A^(x,x')pi(x f )A^(x',x) , (8) 

by means of bilocal averaging operators A% (x, x') and Ap(x',x). The averaging scheme is 
covariant and linear, (ap^ + bq^) = a(pg) +&(q«), a, b G R, by construction and the averaged 
object pp keeps the same tensorial character as pp. 

Let us suppose that the bilocal functions A%(x,x f ) and Ap(x',x) are defined locally on 
an open subset U C A4, x,x' G U. In the following sections their algebraic and differential 
properties are formulated and analyzed to show that averaging operators with such properties 
do exist and also to find out the corresponding properties of the averages 



4 Algebra of the averaging operator 

The following algebraic properties, which are a formalization of the properties of the space- 
time averages in macroscopic electrodynamics using the language of bilocal operators, are 
required to hold p,@,[TJ. 

Property 1 The coincidence limit of A% is 

IrniA^&x') = 5% . (9) 



Property 2 The operator A% is idempotent 

A%(x,x')A^(x',x") = A«„(x,x") . (10) 

These two properties imply that Ap'(x', x) is the inverse operator of Ap,(x, x'), Ap,A^ = 

5™ and Ap,A£ = ftp,, and that the average tensor Pp(x) takes the same value as the original 
tensor Pp(x), Pp{x) = Pp(x), when the integrating region S is chosen infinitesimally small, 
or tends to zero. This implies an additional algebraic property of the averages (||), namely, 
that the averaging procedure commutes with the operation of index contraction. 

The idempotency fliCf ) of the averaging operator A% is designed to provide the idempo- 
tency of the averages (El) in macroscopic gravity 
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Indeed, let us consider for simplicity a vector v a (x) and calculate the twice averaged value 
v a (x) around the same point x G U by using the definition (||) 



T(x) = ^-j (-L / ^(s,sO^,(^s>^V)^^''J ^rfV (12) 



where S x / is an averaging region around the point x' G U. By applying the idempotency 
condition fliCf ) the expression (|i~2"D takes the form 

1 r ( 1 



/fx) 



A%,(x,x")v' y ''(x' , )^j"d n x''\ \pg'd n x' . (13) 



Now, if the term in parentheses were independent of x', we could take it outside the integral 
with respect to the variable x' to get the resulting property (|il"|) . This term in (0), however, 
depends explicitly on x' in the integration region E x /. Similarly as in macroscopic electro- 
dynamics, the idempotency (|ll]) of the averages (||) follows under either of two additional 
assumptions p5| : (nil) the averaging region E x / is the same region E x with the supporting 
point at x' G S x and the average value Pp(x) does not depend on the choice of a supporting 
point within a chosen region; (Hi") the averaging region E x / is a different neighbouring re- 
gion and the average value p%(x) remains the same if evaluated over any neighbouring region 
containing a neighbourhood of x. As it was emphasized in Section 2 this is a fundamen- 
tal property of any physically reasonable classical averaging procedure and such stability of 
averages comes from the basic principles of averaging. 

A microscopic field to be averaged is supposed to have two essentially different variation 
scales^, A and L, satisfying 

A << L , (14) 

and an averaging region must be taken of an intermediate size d such as 

\«d«L (15) 

so that the averaging effectively smooths out all the variations of the microscopic field of the 
scale A. It is implicitly assumed in every averaging scheme (for example, in the averaging 
schemes applied in classical physics in Minkowski space-time, such as in hydrodynamics 
and electrodynamics - see Section 2) that the result of the averaging is insensitive to the 
choice of the supporting point within a fixed averaging region and it is independent of the 



choice of integration (averaging) region itself provided the scale d satisfies the condition ( J15|) . 
This means that the microscopic averaging region is considered as a single "point" for the 
macroscopic averaged field. Such regions have been called "physically infinitesimally small" 



by Lorentz p6 



The trouble is, however, to perform a mathematically rigorous proof of this fact. This 
problem already appears in the usual averaging procedures in flat space-times and, to the 



best of our knowledge, remains unsolved (see, for example, a discussion in [24 for the case of 
hydrodynamics). Furthermore, it is not clear whether the property (|TT|) has to be interpreted 
as an exact one, or it is rather an approximate property. A satisfactory formal analysis of 
the idempotency property of the volume averages is still lacking, though it is extensively 
usedQ This problem certainly deserves further examination [2^ . 

4 The case of two scales is discussed here for the sake of simplicity. Of course, very often there is a 
hierarchy of scales, in which case the arguments are applied for each couple of scales satisfying ( |l4| ) to be 
micro- and macroscoDic ones. resDectivelv. 
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5 Structure of the averaging operator 

In the previous section we discussed two algebraic Properties [I] and |2| of the bilocal operator 
Api(x, x') and the physical motivations which make these two conditions plausible. In this 
section the consequences of requiring @ and fllPl) in the structure of the bilocal operator 
are analyzed. The following main theorem holds. 

Theorem 1 The bilocal operator Ag' (x', x) is idem-potent if and only if it is factorized. That 
is to say, for a bilocal operator A% (V, x) Properties^ and\^ are equivalent to 

A a ;(x',x) = Ff(x')F-^(x) (16) 

where Ff(x) is a set of n linear independent vector fields and F _1 ^(x) is the associated dual 
1-form basis, i = 1, n, and summation over i is carried out. 

Proof. The inverse implication is checked trivially, so let us consider the direct implication 
only. Indeed, the partial derivative of the idempotency relation (fj), A^, (x",x')Ap (x',x) = 

A~p (x",x), with respect to x a ' results in the following expression 



Contracting this expression with Ayr (of, x") and A ,{x,x') and using both (|J) and fllTf ) one 



AZr ta r(x", x')A a p (x', x) + Al, (x", x')Al a r{x' } x) = . 
gets 

Ayrix^X^A^Ax",^) + A^' , x) A^, (x , x') = . (17) 

Now, the first summand does not depend on x and the second one does not depend on x", 
which means that 

A$ t Ax\x)Al l {x : x')=B^, al {x'), (18) 

where B^', a , are arbitrary functions depending only on the variable x'. The relation flllf) 
immediately implies the following identity for the functions B°f, a ,(x'): 

B li'[<r',5'\ + B p'[a' B ^_8'] = > ( 19 ) 

where the square brackets denote, as usual, antisymmetrization and the underlined indices 
are not affected by antisymmetrization. Considering now (|T8| ) as a set of linear partial 
differential equations for the unknowns Ap (x', x) (x being a parameter) 

A$stf,x)=B* a ,tf)Attf,x) , (20) 

one observes that they are always integrable because B^ al (x') satisfy the identity ©. The 
general solution of the system of equations (|20"D is of the form 



A%{x',x) = F*{x')Hi{x) 

where F" (x') 7 i = 1, ...,n, are n 2 independent solutions of (p0|) and H^(x) are n 2 arbitrary 
constants of integration (thereby depending on x). It only remains to impose the coincidence 
limit property @ to set HUx) = F _1 ^(x) and the proof of the theorem is completed. □ 

problems related with idempotency is replaced by the necessity to prove the ergodicity hypothesis which 
states that ensemble and time (or space) averages are equivalent. Both ensemble and volume averagings 
have their own advantages and areas of aDDlicabilitv in describine nhvsical phenomena. It is important to 
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6 Differential properties of bilocal operators 

In this section the differential properties of the bilocal operators are summarized. This 
overview is essential for understanding further results concerning the space-time averaging 
scheme. 

In order to obtain the averaged fields of the geometric objects on Ai one needs to assign 
an averaging region E^ to each point x of U C Ai, where the averaging integral (^j) is to 
be evaluated. Furthermore, to calculate directional, partial and covariant derivatives of the 
averaged fields a law of correspondence between neighbouring averaging regions must be 
defined ||14H,fl6l,[2l]. A natural way to define such a correspondence is to relate averaging 



regions by Lie-dragging, or mapping of a region into another along a vector field (see Section 
2), by means of another bilocal operator Wp (x',x) which is also assumed to satisfy the 
coincidence limit property @. 

To derive the commutation formulae between the averaging and the derivation, one should 
define first the directional derivative of an average field p~g(x) along a vector field £ = d/dX, 

±Fp{x) = Km o ^\p a e(x + Ax)-p^x)] (21) 
where Ax a = £ a (x)AA. Let us define now the shift field for every point 



14 , ie , 2i 



S a '(x',x) = W${x',x)f{x) . (22) 

Now, the averaging region E^+Az associated with the point x + Ax is obtained by Lie- 
dragging the averaging region a parametric length A A (the same for all x' G E^) along 
the integral lines of the field S a . 

The coordination bivector Wp allows the construction of the shift vector S a for any 
averaging region and any vector £ Q (the shift vector at the supporting point for that region). 
By choosing n such linearly independent vector fields £f and shifting averaging regions along 
them one can build a covering of the manifold with one averaging region associated to every 
x G 1A C Ai. This procedure is a formalization of the condition (i) of Section 2. 

As a consequence of the definitions (|6]), ( 21) and (p2j ) one can obtain the following formula 



for the commutation of partial differentiation and averaging jnj . |Tj| : 

n,x = (A^^Wi) + (p^VVL) - - (3aP» + (p a A) . (23) 

Here, S£ x are the so-called structural functions, Sp y = A"(Ap n + A e p a ,W° ) and W|. e / is 



the divergence of the coordination bivector, W|'. e / = Wjf e ' + (In yj— g)^W(. In an affine 
connection space with connection coefficients (a Riemannian space is considered in ||21|| ) 
the same formula (|23|) remains valid for covariant derivatives, where partial differentiation is 
replaced by the covariant one and the divergence W e x . e , is replaced by Wf. e , = W{ e , + IV W{ 
(with r e / = r"' a ,). The commutation formula ( p3| ) has a very transparent meaning: the first 
term in the right-hand side is the average value of the derivative of weighted by Wp , 
the last two terms are due to the non-triviality of the averaging operator Ap\ while the 
second and third terms describe the effect of a non-trivial averaging measure in (|6|) and the 
variation in the value of the volume for different regions, respectively. Indeed, the change in 
the volume Vs of an (averaging) region E C A4 Lie-dragged along a vector field £ is given 
by (see, for example, |^0fl,[P31) 
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which can be written in terms of partial derivatives after using ( p2|) as 

Vjw = (W# a ,)V s . (25) 

The expression (|23| ) is the most general version of formulae (0) for arbitrary averaging and 
coordination operators. 

A fundamental problem one has to face now is whether it is possible or not to define a 
unique covering of the manifold A4, like in Minkowski space-time (see Section 2). In other 
words, given a microscopic tensor field Pp(x), the goal is to determine a uniquely defined 
averaged tensor field p^(x) with reasonable analytical properties in its dependence on the 
supporting point. Furthermore, we must deal with the problem that the averages (|>D de- 
pend functionally on the averaging regions Pp(x) = p%(x)[Yi\ and applying the commutation 
formula (|23|) to find averaged equations would, in general, bring volume dependent terms 
into them, which is undesirable. In the case of Minkowski space-time the uniqueness of the 
averaged field with the proper analytical behaviour @ and the parametric dependence on 
the averaging region volume are ensured by the condition (ii) (see Section 2) which sets the 
averaging regions to be typical, of the same shape and volume, Vg = const and At = const, 
and shifted along the Cartesian coordinate lines. To arrange similar properties for the gen- 
eralized averages (^) one must look for some specific conditions on the coordination operator 
. The following remarkable theorem holds |ljj],|l6| vers i° n 101 Riemannian manifolds 



has been given in pTf ). 



Theorem 2 In the averaging region coordination by the bivector VVg' on an arbitrary dif- 
ferentiable manifold, it is necessary and sufficient to require 

Wfp,H + = , (26) 

for the average tensor field p2(x) to be a single valued local function of the supporting point 
x onU C M, 

Pp,M = • ( 27 ) 
Proof. The sufficiency of ([26]) is easily proved by calculating the antisymmetrized second 
partial derivative of (^) which, after using (p6|), gives (27]). The necessity follows from the 



analysis of the condition (|27|) after taking into account that the averaging regions and the 
tensor field p^ are arbitrary. □ 

Geometrically, (^) means that, given an averaging region E, the region S obtained 
by transporting S along an infinitesimal parallelogram constructed from two commuting 
vector fields £ and ( according to the law ( p2|) coincides with the original region, E = E. 
This is a highly non-trivial property which allows to construct a covering of the manifold 
with an averaging region attached to every point in the manifold, thus generalizing the 
corresponding part of the condition (ii) in Section 2. In the formalism of bilocal exterior 
calculus the condition ( PB[ ) reads that the operator is biholonomic, which means that 
the bilocal coordinate 1-form basis W a ' = Wp'dx 13 has vanishing biholonomicity coefficients 
in the bilocal Maurer-Cartan equations |16[ . 

It should be noted here that formula (^) is analogous to formula (§) in both form and 
meaning. 

Another condition on the coordination bivector is the requirement that the Lie- 



dragging of a region is a volume-preserving diffeomorphism |l4|,|[L6| ([21] for Riemannian 
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which means that the averaging regions do not change the value of the volume when shifted 
(coordinated) along a chosen vector field £ according to (|22|) . This generalizes the corre- 
sponding part of the condition (ii) in Section 2. 

Thus, the condition fl26|) states that the average tensor field is a single valued local 
function of the supporting point x (ffil )- Adding the condition (^Sf), the average tensor field 
does not depend explicitly on the value of the region volume V = Vs, and V itself is a free 
parameter of the theory. Given a microscopic tensor field Pp(x) on Ai, the average tensor 
field pp (x) is therefore uniquely defined on U C .M and can be handled within the framework 
of standard differential and integral calculus. Requiring additionally that the two bivectors 
Ap and Wp coincide 

A* = Wg (29) 

the first term in the commutation formula (|23"D becomes exactly the average derivative. Using 
all conditions (fffi), fl2"8"|) and (P9"[), the commutation formula acquires a remarkable simply 
form []14|,|I6| 

^ i7 = (p^ + p^,>V 7 a '). (30) 

The corresponding analogues of this expression for covariant differentiation are obtained 
by replacing partial derivatives by covariant ones. To obtain the expression for directional 
derivatives we must contract this expression with a vector £ 7 and insert the vector field 
S a from fl22"|) in the second term of the right-hand side of ([3(]). Formula fl3"0|) generalizes 
formulae (0) and it can be easily shown [16[] to become exactly (§) if Wf' = 5£ and the 
volume n-form e is standard, that is (In s/—g) >e / = and e = dx 1 A ... A dx n P^] , [ fZ3]] . In 
Section 8 one can find more details on this particular case, see formulae (f46l) and (B3|). 



7 Existence theorems 

The differential conditions (|26D and ( p8|) together with the algebraic conditions (§), ([10|) 
and ( |29| ) are to be considered as a set of partial differential and algebraic equations for the 
unknown functions Wp . Provided a solution for the system is found, the existence of such 
operators Ap' and Wp (and therefore of the averages with the above described properties) 
is proved. Theorem [3] has revealed the structure of the operator Wp obeying the algebraic 
properties (|9|), ([[(]). The following theorem gives the general solution of ( p6|) (in fl4|| ,[|T6|j 



and |2T| for Riemannian manifolds, the same theorem has proved a solution of fl26|) with (29 ) 



sought in a factorized form (]T5|) , and now, with Theorem [I] taken into account, it gives the 
general solution). 



Theorem 3 In an arbitrary n-dimensional differentiable manifold the general solution of 
the equations 

'V 7 ] T ¥V [P,S' i V ' \ 
for idempotent bilocals Wp'(x',x) is given by 



< 1 + WSV>Vf; = , (31) 



W a p'(x',x)=f?(x')f- 1 p(x) (32) 
where f° L {x)d a = fi is any vector basis satisfying the commutation relations 
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with constant structure functions (anholonomicity coefficients) C^, 

C% = const. (34) 
Proof. Due to Theorem [l] the idempotent bivector W% (x', x) has the factorized form 

m 

W$\x\x) = F?{x')F-%x) . 
Inserting this expression into ([H]) one finds 

K \x')F-\^x) + F«Ux')F-^(x)F?(x')F-^(x) = . (35) 
Now, using the expression relating derivatives of a vector basis and its dual 1-form basis 



in equation fl35]), it becomes 



= -FfF-^F-* (36) 



2F[F° k F-\ a] (x) = 2Ff F£ F' 1 ^ (x') = C\ k , (37) 

which is exactly the expression for the anholonomicity coefficients C\ k in terms of a 1-form 
basis. Equations fl37[) say that any vector basis F"(x) = f^{x) such that its corresponding 
dual 1-form basis f~ 1% = f~ 1% a dx a satisfies the Maurer-Cartan equation 



1 



dT 1 ' = -~CLf- lj A f~ lk (38) 



3 



with constant anholonomicity coefficients C\ k = const, is a solution of (|3i"l). The class of 
all such bases constitutes the general solution of equations ([[!]) and therefore the general 
solution for the coordination bivector satisfying ( pB"l) and (|29|) . □ 

The next theorem proves the existence of solutions for the equation (|28|) within the class 
of bivectors satisfying (^) p5f (its version for a particular subclass of (j32|) , see Section 8 
and Proposition 1 below, has been given in |I3]],[[nj, and |^T| for the case of Riemannian 
manifolds) . 

Theorem 4 In an arbitrary n- dimensional differentiable metric manifold (M, g a p) with a 
volume n-form there always exist locally volume-preserving bivectors Wg (x',x) of the form 
(HU with (j§£) satisfying ([2$). 

Proof. For a bivector Wp'(x',x) = ff {x')f~ l% p{x) equations (p8[) read 

fZa> = f&> + (^V=g),«>f? = o. (39) 

Writing (|3|) as 

-mv + On^-O, (40) 



and using the definition of the anholonomicity coefficients in terms of 1-form basis (see (|3 
above) this equation can be rewritten as: 
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Now it is sufficient to show the integrability of this quasi-linear partial differential equation 



to prove the theorem. Taking into account (|34]) and the property C^Cjy = which holds in 
that case, one can easily show that ([IT]) is always integrable on an arbitrary n-dimensional 
differentiable manifold with a volume n-form. □ 

It should be stressed here that the conditions Ck = const due to Theorem ^| are essential 
for the proof of Theorem |] and they guarantee the local existence of n linear independent 
divergence free vectors, the result holding for both orientable and non-orient able manifolds 
(by using the so-called odd volume n-form ||28|| , |f29H )p|. 

Another important point is that n vector fields /" satisfying (|33|) and (]34]) define a 
finite dimensional Lie group on the averaged manifolds which is directly related with the 



symmetries of such manifolds [[31 



These two theorems prove the existence of solutions for the set of equations ([|), fllPp, 
fl26|), fl2~8|), and (]29|), and therefore prove the existence of the bilocal operator (x', x) with 
the corresponding algebraic and differential properties (^) and fl3Tj| ) for the averages (j^). 

Now a particular subclass of the operators ([32]) with fl3~j| ) will be considered to analyze 
some additional properties and to reveal the functional structure of the subclass. 



8 The proper systems of coordinates 

As it has been emphasized in the Introduction and Section 2, the space-time averaging 
procedure adopted in electrodynamics is essentially formulated in Cartesian coordinates and 
all its properties are shown by exploiting the exceptional character of this coordinates (see 
Section 2 for details). The covariant formalism developed for the averages (J6|) in macroscopic 
gravity (Sections 3-6) generalizes the averaging scheme of macroscopic electrodynamics for 
arbitrary n-dimensional differentiable manifolds and while keeping covariant properties which 
are analogous to those in electrodynamics. 

Let us now consider the macroscopic gravity averaging scheme for a particular subclass 
of operators fl32|) with (34). This particular subclass admits a special coordinate system in 



which the averages and their properties have especially simple form and meaning. Such a 
coordinate system is an analogue for macroscopic gravity of the Cartesian coordinates in 
Minkowski space-time. 

Let us hereby restrict the class of solutions of the equations ([H]) to the subclass satisfying 

[fii fj] = , (42) 

that is C*j = 0. In this case the vector fields jf constitute a coordinate system and there 
always exist n functionally independent scalar functions <fi l (x) such that the vector and 
corresponding dual 1-form bases are of the form 

RWfi)) = %, r 1< M(x»)) = %- (43) 



n linear independent divergence free vector fields on an n-dimensional differentiable manifold with a 
volume n-form with structure functions (|33| ) can be shown to exist locally iff C\ k i = 0. This condition is 
fulfilled also globally on parallelizable manifolds (both orientable and non-orientable) |j(J. This shows that 
Theorem [| is valid globally for such manifolds (a manifold is called parallelizable if its tangent bundle is 
trivial") . 
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Thus, the bilocal operator W« (x',x) becomes 

Being functionally independent, the set of n functions cj) l {x) can be taken as a system of 
local coordinates on the manifold M. ]14[ , [16|] (|2l|] for Riemannian manifolds). 

Definition 3 A coordinate system {</>*} defined by n scalar functions (ff = <ft l (x) in jff^j ) will 
be called a proper coordinate system. 

The usefulness of this definition is motivated by the fact that in a proper coordinate 
system the bilocal operator Wp (x', x) takes the simplest possible form 

W;(0',0) = = < = S), (45) 

where the bilocal Kronecker symbol dp is defined as Sp* = SfSp. The definition of the 
average (^) acquires a remarkable simple form (closely resembling the space-time averages 
of macroscopic electrodynamics ([!])) when written using a proper coordinate system 

P}(0) = TT f P}(0OV-^W • ( 46 ) 

Theorem || has proved the existence of volume-preserving bilocal operators fl32|) and fl34|) , 
and therefore the existence of solutions of (28) for operators Wp of the form (f44[) from this 
theorem. It is useful, however, to prove it here independently due to the above mentioned 
importance of the proper system of coordinates for the macroscopic gravity averaging scheme. 
The following statement holds in this case |TTJ|,|T^] (|f2T| for Riemannian manifolds). 

Proposition 1 In an arbitrary n-dimensional differentiable metric manifold (Ai, g a p) with 
a volume n-form, there always exist locally a set of n scalar functionally independent func- 
tions <p l (x) such that the corresponding coordination bivector (jQ) satisfies condition (jlffi). 

Proof. Let us write down equations (|28| ) for the bivector Wp' ([44]) 

% + (\nJ^)\ ^=0 (47) 



d(p l d(j>> dx c 

in terms of the unknowns <fi l {x). These equations are equivalent to the system of quasi-linear 
partial differential equations (compare with (HTJ) ) for f~ 1% a (S3) 



- f?'(x'(<P k ))r%(<f> k ) + (hV-sW))) . = (48) 

which are always locally integrable on an arbitrary n-dimensional differentiable manifold 
with a volume n-form. Then, in accordance with (f£J) <fi l (x) are solutions of the equations 
<f) l a = f~ ll a , which are also always locally integrable because the integrability conditions 
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Corollary 1 Any proper coordinate coordinate system such that the corresponding bivector 



satisfies condition ( \2q ) is a volume-preserving system of coordinates, 



lnJ-g{<f) k ) = , or, g(<p k ) = const . (49) 



It should be noted here that all arguments concerning non-orientable manifolds and global 
existence given after Theorem |] apply here as well (see Section 7) . 

For the case of (pseudo)-Riemannian manifolds, Corollary [l] states that in a proper coor- 
dinate system the Christoffel symbols (which are T^ a = (ln-^/ = g) i( g due to equi-affinity 
of Riemannian manifolds) vanish 



T i (<j>) = Ti j = [\nJ-gm] =0. (50) 



Another useful characterization of the volume-preserving coordinates, in addition to (f49|) , 
or (|50"D , can be obtained in terms of the expansion of the vector fields tangent to the coor- 
dinates lines. Defining a vector X(i) tangent to a coordinate line 0* as 

*« = 4>^ = 44p ( 51 ) 

it is immediate to find that in the proper coordinate system {</>*} the expansion divx(i) of 
the vector field ( 51 ) is 



4):j = [H-^)) . = , (52) 

so that the condition of vanishing expansion for the tangent vector fields X(i) i s equivalent 
to the definition of volume-preserving coordinates (|49"|). 

In addition to very simple and transparent forms of the coordination bivector fl4"4] ) and 
the averages (fR]), the volume-preserving proper coordinate systems also allow a remarkably 
simple expression for the commutation between partial differentiation and averaging. Indeed, 
expression ( |30"D becomes in the proper coordinate system 

9 :m) = (^P i 0)) , (53) 



which is exactly the same commutation formula for as in the averaging scheme in Minkowski 
manifolds (@). 

Let us now study the functional structure of the class of the volume-preserving coordi- 
nates to understand how large it is and how much freedom for coordinate transformations 
it contains. The following Proposition reveals the structure of the class. 

Proposition 2 The class of volume-preserving coordinate transformations on an arbitrary 
n-dimensional differentiate metric manifold (Ai, gap) with a volume n-form, contains (n— 1) 
arbitrary functions of n arguments and one arbitrary function of (n — 1) arguments. 

Proof. Due to Proposition [I] the volume-preserving coordinates are characterized by the 
condition (f49D - The class of coordinate transformations y a = y a (x tl ) which preserve this 
condition, 
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is denned by the system of partial differential equations |3 

d 2 y a dx^ 



dx v dx^ dy c 



. (55) 



Equations ( |5"5"D can be easily written as equations for the Jacobian of the coordinate trans- 
formation y a = y a (x' 1 ), 



_d_ 

dx h 



det 



dy a 
dx^ 



, (56) 



which have the general solution 



**(g)=C (57) 



where C is a non-vanishing constant to ensure invertibility of the coordinate transformation. 
Expanding the Jacobian ([57]) by its first row we get 

|; ? '(y 2 ,...,,») + ... + |lfV,...,y") = C (58) 

where £, u {y 2 , y n ) are the corresponding minors of the determinant, depending therefore 
only on (n — 1) functions y 2 , ■ ■■,y n - These (n — 1) functions are arbitrary functionally inde- 
pendent functions of x^, and can be taken as explicit functions of x M , = £ y (x M ). Then 
equation fl5"8|) reads 

This is a linear partial differential equation for the unknown y l (x^). Its general solution is 
known to be 

y V) = Vi(*") + *i(*V), k n {x»)) (60) 

where yp{x^) is a particular solution of (|59"D , k (x**), A = 2,...,n, are n — 1 functionally 
independent first integrals satisfying 

r^ = o, (6i) 

and h is an arbitrary function of (n — 1) arguments fc" 4 . Thus, the general solution of the 
equations (|55|) defining the class of volume-preserving coordinates fl54|) contains (n — 1) 
arbitrary functions y A (x fi ) of n arguments x M and one arbitrary function h(k A ) of (n — 1) 
variables fc A (x M ). □ 

The set of proper coordinate systems forms quite a big class and this functional freedom 
may be used to specify additional properties of the averages or (p[6|), when necessary. 

Choosing different proper coordinate systems (ft will give different average fields fl46f) 
of a given microscopic tensor field p^(x). In general, the averages p*(0) and Pj(<ft) calcu- 
lated in the proper coordinate systems (ft and (ft are not related by a tensorial law under 
the transformation (f> % = <f) l (<ft), nor are the operators Wj(0', (ft) and Wj (</>', 0) related by 
a tensorial transformation with each of them being <5*- (^) in its own proper coordinates 
system. It should be noted here that averages @, and (|46l), are obviously tensorial with 
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fl32]), or (|44|), of the bilocal operator Wp'(x',x) itself, which involves a functional freedom 
in changing the functions f"{x), or <p l (x). This "non-tensorial" property is very natural, 
indeed, for it states the exceptional character of the proper coordinate systems for obtaining 
the simplest and most transparent form of the averages and the averaging and coordination 
operators. It closely resembles, on the other hand, the definition of averages in macroscopic 
electrodynamics as it was noted above, and the exceptional character of Cartesian coordi- 
nates used in that averaging procedure. The class of the proper coordinate systems on an 
arbitrary differentiable metric manifold is a natural counterpart of the Cartesian coordinate 
system on a Minkowski manifold. The property they share in common is that both are 
volume-preserving. 

There is, however, a special subclass within the class of volume-preserving coordinate 
transformations described in Proposition [2] which keeps the bilocal operator (^) and the 
averages (f46|) covariant. 



Proposition 3 The class of transformations (f)' 1 — » (f) 1 which keeps the bivector W^' and the 
averages (||) covariant within the class of proper system of coordinates is 

ft = A}0 J + a} (62) 

where A* and a 1 are constant. 

Proof. The proof is straightforward. Consider two proper coordinate systems <p l and 0*, 
and require the bivectors Wj(0',0) and Wj(0',0), each of the form (pl|) in its own proper 
coordinates, to be equal 

dx a ' dcj) 1 dx a ' dft 

Since the coordinate system x a in equation fl63| ) is proper by assumption it must be either 
<p % or cf) 1 . Choosing <p l one immediately obtains 

gi = dft_ dft_ 
3 d&tfwWM*)' 



This last equation has as general solution (p^) . It is easy to show that the class of transfor- 
mations between proper coordinates keeps the averages ( pf) also covariant. □ 

Due to Proposition |, if the manifold (A4, g a p) is chosen to be a (pseudo)-Riemannian 
spacetime, the averages ( ^6| ) defined in proper coordinates are Lorentz tensors exactly like 
the averages in Minkowski space-time (see Section 2). 



9 Conclusion 



Thus, the covariant averaging procedure for objects (tensors, geometric objects, etc. 
framework of macroscopic gravity 



•0,0 



in the 

is a natural generalization of the space-time 
averaging procedure of macroscopic electrodynamics. It gives a covariant formulation of the 
conditions on the averages, which provide them with natural algebraic and analytical prop- 
erties. A wide class of averaging and coordination bilocal operators satisfying all the proper- 



9 CONCLUSION 



17 



spaces. The class of proper coordinate systems, analogous to Cartesian coordinate system of 
Minkowski space-time and generalizing them, gives the simplest and most transparent form 
of the averages. The averaging procedure as it is formulated allows a large functional free- 
dom which is incorporated in an elegant way and can be used to arrange additional specific 
conditions for the averages. 
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